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ABSTRACT 

In the first paper of this series we use the publicly available code Gim2D to model 
the r- and i-band images of all galaxies in a magnitude-limited sample of roughly 
1800 morphologically classified galaxies taken from the Sloan Digital Sky Survey. The 
model is a concentric superposition of two components, each with elliptical isophotes 
with constant flattening and position angle. The disk luminosity profile is assumed 
exponential, while the bulge is assumed to have a de Vaucouleurs or a Sersic profile. 
We find that the parameters returned by Gim2D depend little on the waveband 
or bulge profile used; their formal uncertainties are usually small. Nevertheless, for 
bright galaxies the measured distribution of b/a, the apparent disk flattening, deviates 
strongly from the expected uniform distribution, showing that the "disk" identified 
by the code frequently corresponds to an intrinsically 3-dimensional structure rather 
than to a true thin disk. We correct approximately for this systematic problem using 
the observed statistics of the b/a distribution and estimate, as a function of absolute 
magnitude, the mean fractions of galaxy light in disks and in "pure bulge" systems 
(those with no detectable disk). For the brightest galaxies the disk light fraction is 
about 10% and about 80% are "pure bulge" systems. For faint galaxies most of the 
light is in disks and we do not detect a "pure bulge" population. Averaging over the 
galaxy population as a whole, we find that 54 ± 2% of the local cosmic luminosity 
density at both r and i comes from disks and 32 ±2% from "pure bulge" systems. The 
remaining 14 ± 2% comes from bulges in galaxies with detectable disks. 

Key words: Survey - galaxies: photometry - galaxies: bulge-to-disk decomposition 
- galaxies: fundamental parameters - galaxies: spirals - galaxies: ellipticals - galaxies: 
bulge - galaxies: disk 



1 INTRODUCTION 

Observations show that galaxies typically have two compo- 
nents with diflferent photometric and dynamical properties: 
a thin, rotationally supported stellar disk which often also 
contains gas and dust, and a bulge, a spheroidal or ellipsoidal 
component made purely of stars. Most galaxies possess both 
structures, but some, the ellipticals, have no significant disk, 
while others, late-type spirals and irregulars, have little or 
no bulge. Since these two components presumably had dif- 
ferent formation paths, their relative importance must be a 
fundamental clue to how each galaxy formed. Among the 
disk galaxies themselves there exists a dichotomy between 
those rich in gas, dust and the accompanying star forma- 
tion, spiral and irregular galaxies, and those in which this 
activity is virtually absent, lenticular galaxies. 

It is known that galaxies of differing morphology are 
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segregated according to environmental density IIDressleii 
1980) with the low-density field composed largely of spirals 
and irregulars and the densest regions of clusters composed 
of lenticulars and ellipticals. Many interpretations have been 
proposed for this observation, ranging from initial condi- 
tion biases which imprint the differences at birth, to gravi- 
tational, gas dynamical or radiative processes through which 
galactic environment affects later evolution. Recent studies 
using the large, complete samples provided by the Sloan Dig- 
ital Sky Survey (SDSS, see below) show the distribution of 
galaxy mass to depend quite strongly on environment, ap- 
parently reflecting an initial condition bias, while for galax- 
ies of given mass, star formation properties depend much 
more strongly on environment than do structural properties, 
suggesting that external processes primarily affect the ga s 
component from which stars form jKauffmann et a'l]|2004^ . 

Understanding the chronology of bulge and disk forma- 
tion by analysing the relative contributions of these two com- 
ponents at different cosmological epochs is a fundamental 
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goal of observational cosmology. Such data provide impor- 
tant constraints on competing scenarios of galaxy formation 
and evolution. In models where bulges form first and disks 
are added later no close correlation is expected between the 
two formation phases; the ratio of spheroid luminosity to to- 
tal luminosity measures the the efficiency of the first burst 
of star formation relative to later slow accretion (Fall & Ef- 
stathiou 1980; Kauffmann, White & Guiderdoni 1993). If 
instead bulges form out of disks by secular evolution, then 
a stronger correlation between their properties may be ex- 
pected (Norman et al.l996; Bournaud & Combes 2004). 

The relative contributions of bulges and disks to the 
stellar content of galaxies of differing mass and also to the 
overall stellar content of the Universe are clearly important 
quantities which should be reproduced within any viable 
picture of galaxy formation. It has long been known that 
the massive galaxy population is dominated by light from 
spheroids, w hile disks tend to domin ate in lower mass pop- 
ulations (e.g. ^Efs tathiou et al.l il982h '). but there have been 
no recent quantitative studies for the galaxy population as a 
whole. The most commonly cited estimates of the fractions 
of all stars in bulges and disks trace back to the work of 
ISimien fc de Vaucouleur j lll98d) and ISchechter fc DressleJ 
These results suffer from poor statistics and are 
based largely on photographic material and on visual "de- 
composition" of the light of each galaxy into its bulge 
and disk fractions. With large photometric surveys, linear 
CCD detectors and quantitative decomposition techniques 
it should clearly be possible to do much better. This is the 
goal of the present paper and subsequent papers in this se- 
ries. 

We apply a modern two-dimensional morphological de- 
composition algorithm to the images of a magnitude-limited 
sample of relatively nearby and bright galaxies with photom- 
etry available from SDSS. In this paper we use the results to 
study the overall contribution of disks and bulges to the light 
of the nearby galaxy population both as a function of intrin- 
sic galaxy luminosity and for the population as a whole. In 
section 121 we describe the selection of the sample we use. In 
section 13 we describe the decomposition algorithm and the 
parametric functions it uses to fit the galaxy images. We 
also describe how the data were processed. In section 3] re- 
sults are compared for different photometric bands and for 
different parametric fitting functions in order to explore the 
robustness of the results. Deviations of the images from the 
best fit models are also quantified. In section|21we study the 
distribution of derived parameters for those galaxies which 
could be successfully fitted. We also identify and correct for 
a serious (and previously known) systematic which results in 
the assignment of significant disks to bright galaxies where 
disks are in most cases actually absent. This section presents 
our principal results, estimates of the fractions of the total 
light of all galaxies of given luminosity which are in disks, 
in bulges, in galaxies with no detectable disk or in galaxies 
with no detectable bulge. Combining these with previous 
measurements of the galaxy luminosity function allows us 
to obtain corresponding fractions for the galaxy population 
as a whole. A final section discusses and summarises these 
results. 

Throughout this paper, unless otherwise stated, we as- 
sume a Hubble constant of Ho = 70 km Mpc~^ and an 
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Figure 1. The distributions of galaxies in morphological type 
(upper panel) and in extinction-corrected Petrosian apparent 
magnitude in the r band (second panel) are shown for the 1862 
galaxies in our sample. The redshift (third panel) and correspond- 
ing absolute magnitude (lower panel) distributions are given for 
the 1550 galaxies in the sample for which we have spectroscopic 
information. 



Qm = 0.3, SIa = 0.7 cosmology in calculating distances and 
luminosities. 



2 OBSERVATIONAL DATA 
2.1 Galaxy sample 

In June 2 01 th e Sloan Digital Sky Survey (SDSS; 
I York et alJ (l2 000l)) released its Early Data Release 
(EDR lStouehton et all ^00 f)), roughly 462 square degrees 



of imaging data collected in drift scan mode. The imag- 
ing is conducte d in t he w— , q— , r— , i— and 2— bands 



Gunn et alJ 199^: iHogg et alJ 



ifFukugita et al.' h99&\ : 

(2001); Smith ct al. ( 2002 h IPier et al.1 j2003() l. The reader 
is referred to llvezic et al ^QQ^ioi details on the phto- 
metric quality assessment. The EDR contains around a mil- 
lion galaxies distributed within a narrow strip of 2.5 degrees 
across the equator. As the strip crosses the galactic plane, 
the data are divided into two separate sets in the North 
and South Galactic caps. The SDSS has the ambitious goal 
to image a quarter of the Celestial Sphere at high Galac- 
tic latitude as well as to obtain spectra uniformly for all 
the brighter galaxies. For the present project this has the 
advantage, in comparison to previous work, of having uni- 
form photometry and spectroscopy over a much larger area, 
permitting a major improvement in sample size and homo- 
geneity. 

In the following analysis we are using a sample of 
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Subsamples 




Sample 


< T < 1 
E 


1 < T < 2 
SO 


2 < T < 3 
Sa 


3 < T < 4 
Sb 


4 < T < 5 
Sc 


5 < T < 6 
Sd 


T = 6 
Irr 


T = -1 
unclassified 


Total 


Photometric 


487 


417 


313 


312 


232 


48 


25 


28 


1862 


Spectroscopic 


413 


363 


272 


262 


197 


42 


16 


23 


1588 



Table 1. Visual classification of our photometric and spectroscopic samples into morphological subsamples. 



galaxies defined by the Japanese Participation Group (JPG, 
iNalcamura et alJ |^03)). This is a homogeneous sample ob- 
tained from the northern equatorial stripes of the SDSS 
EDR. The region of the sky covered is an area of 229.7 square 
degrees at 145.15° < a < 235.97° and -1.27° <S < 1.27°. 
The sample is limited to bright galaxies with r < 15.9 after 
Galactic reddening correction. Eye classifications cannot be 
made confidently beyond this magnitude, and for our cur- 
rent purpose this bright limit has the advantage that the 
galaxies are all large compared to the SDSS point-spread 
function. 

All the 1862 galaxies in the sample were cla ssified by 
eye o n the system of the Hubble Atlas of Galaxies dSandaed 
llQGlT) by JPG scientists using the g band image of each 
galaxy. For each galaxy the final quoted type is the mean of 
4 independent classifications by different scientists. The rms 
of these 4 classifications is also given; they typically agree 
within AT < 1.5. The corresponding numerical classifica- 
tion as d efined in the Third Referen ce Catalogue of Bright 
Galaxies ide Vaucouleurs et al J [199111 is also reported. The 
seven resulting subsamples separate galaxies according to 
morphology going from ellipticals (E) through lenticulars 
(SO) and early-type spirals (Sa, Sb), to late-type spirals (Sc, 
Sd) and irregulars (Irr). For 1588 galaxies out of our sam- 
ple of 1862 we have spectroscopic information. A summary 
of the distribution of our galaxies across the morphologi- 
cal cl asses is given in Tabled We refer to jNakamura et all 
l2003l) for further details. Further work ex ploring the proper - 
ties of this sample can be f ound in jNakMnuraet all i2004l) . 
iKellv fc McKavl (|20oi) and lFukueita et al.l j20oJP 



2.2 Photometric and spectroscopic data 

Two important quantities used in this paper are taken di- 
rectly from the SDSS database: the redshift and the Pet- 
rosian magnitude. The first is obtained by the spectroscopic 
pipelines idllspec2d (written by D.Schlegel & S. Buries) and 
spectrold (written by M. SubbaRao, M. Bernardi and J. 
Frieman) . A description of the tiling algorithm used to assign 
targets to each pointing is given in Blanton ct al. (2003). 
The SDSS spectroscopic gala xy samples consist of all galax- 
ies brighter than r = 17.77 listrauss et al.ll2002h and of a 
sample of luminous red galaxier TEisengteine^il.ll200 iTl ex- 
tending at r < 19.2. The distribution of galaxies with respect 
to z for our sample is shown in the third panel of Figure 
The second quantity is obt ained by the Photo pipeline (see 
iLupton et alJ (1200 ll . |20o3)) and it is based on a modified 
form of the Petrosian system for galaxy photometry which is 
designed to measure a constant (and large) fraction of the to- 
tal light of a galaxy independent of its characteristic surface 
brightness. Three related quantities also used in this paper 
are Rsq and Rgo, defined as the radii which include respec- 
tivel y 50 and 90 perc e nt of the Petrosian flux in the r band 



c = Rgo/Rso- The second panel of Figure^shows the distri- 
bution of galaxies with respect to their Petrosian magnitude 
after correction for foreground Ga l actic e xtinction using the 
reddening map of ISchleeel et al.l (Il99d) . Such extinction- 
corrected Petrosian magnitudes are used throughout this 
paper. 



3 IMAGE ANALYSIS 
3.1 The fitting algorithm 

We exam ine the structural p roperties of our galaxies using 
Gim2D (Sim ard et al .120021) . a two-dimensional photomet- 
ric decomposition algorithm which fits each image to a su- 
perposition of an elliptical component with a Sersic profile, 
representing the bulge, and a concentric elliptical compo- 
nent with an exponential profile, representing the disk. It is 
important to recognise that this separation into bulge and 
disk is based only on the observed image, and may not corre- 
spond to the "best" decomposition if additional information, 
for example from kinematics of the stars or the gas, is taken 
into account. 

It is well known that the bulge component of most 
galaxies can be well represented by a surface brightness pro- 
file of Sersic form; 



E(r) =Ee ■exp{-6[(r/re) 



l/r. 



1]} 



(1) 



where S(r) is the surface brightness a distance r from the 
centre along the semi-major axis and Ee is its characteristic 
value, the effective surface brightness, defined as the value 
at the effective radius r^. The parameter b is related to the 
Sersic index n and is set equal to 1.9992n — 0.3271 so that r^ 
is the projec ted ra dius enclosing half of the total light iSersid 
(1968; Ciot'til ll99J) . Many authors fit bulges and ellipticals 
with a more specific function, the de Vaucouleurs r^^* law, 
which is obtained by setting n = 4. When fitting Sersic pro- 
files in the following we will assume 0.2 < n < 4. This choice 
is driven by the knowledge that fits to low-luminosity ellip- 
ticals and bulges generally give n significantly smaller than 
4, the standard value for bright ellipticals (de Jong 1996; 
Caon et al.2005). Values of n in excess of 4 are sometimes 
found for cD galaxies but our sample does not contain such 
exceptionally luminous systems. 

Disks are generally well described by an exponential 
profile (corresponding to n = 1) , although non-axisymmetric 
features due to bars, spiral arms and dust lanes can be large. 
We use the standard parametrisation: 



E(r) = Eo • exp(-r//i) 



(2) 



fsee IStoughton et alJ (120021) '). and the concentration index 



where Eo is the central surface brightness and h the disk 
scalelength. 

These laws are purely empirical fitting functions with no 




Figure 2. Examples of the science image, the mask, the model and the residual images for our two fits (de Vaucouleurs plus exponential 
or Scrsic plus exponential) in the r band (from left to right) for galaxies from our seven morphological classes: E, SO, Sa, Sb, Sc, Sd, Irr 
(from top to bottom). 
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strong theoretical justification. Other functions might pro- 
vide equally good fits to a galaxy profile leading to a dif- 
ferent parametrisation for the bulge-to-disk decomposition, 
and possibly different values for global parameters such as 
the bulge-to-total luminosity ratio B/T. 

With this model there is a maximum of twelve parame- 
ters which are adjusted in fitting the galaxy image and that 
we retrieve as output from our decomposition: the total flux 
of the object; the bulge-to-total light ratio B/T, defined as 
the fraction of the total flux in the bulge component so that 
B/T = 1 corresponds to a pure bulge and B/T = to a 
pure disk; the bulge effective radius r^; the disk scalelength 
h; the disk inclination angle i defined so that i — for face- 
on disks and i = 90 for edge-on ones - the disk axial ratio 
is then {b/a)disk = cos(i); the bulge ellipticity e given in 
terms of the bulge axial ratio by e = 1 — ip/a); the bulge 
and disk position angles (hereafter PA) measured clockwise 
from north and allowed to be different; the Sersic index n 
which we sometimes fix at n = 4; the x-y pixel shifts dx and 
dy of the galaxy centre position in the model and science 
thumbnail images; and the background intensity level. 

Additional parameters could be introduced to model 
other features (e.g bars, spiral arms, etc.) but such decom- 
positions become somewhat arbitrary and may not converge 
to unique solutions. Our current choice is standard and we 
found it to be a good compromise between stability of results 
and fiexibility of representation. Notice that for an axisym- 
metric galaxy the position angles of the bulge and the disk 
would be the same. Allowing them to differ makes it possi- 
ble for the code to detect triax;ial bulges or bars in suitably 
oriented galaxies. 

3.2 Image reduction 

We perform our analysis starting from corrected frames of 
area 13'. 52 x 8'. 98 taken directly from the SDSS archive. 
In these large-scale images flat-field, bias, cosmic-ray, and 
pixel-defect corrections have already been applied. The 
pixel-size in these images is 0,396 arcsec. To proceed with 
our fitting, we begin with a list of source positions and apply 
the SExtractor galaxy photometry package version 2.2.2 
jBertin fc Arnout jll99g) to each field to estimate the local 
sky background level at each point and to define the isopho- 
tal area where each object is above the detection threshold 
(we choose a threshold which is higher than background by 
1.5 times the background noise). When SExtractor performs 
galaxy photometry, it constructs a segmentation (or mask) 
image in which pixels belonging to the same object all have 
the same value and sky background pixels are flagged by 
zeros. Our 2-D image fit is carried out on all pixels be- 
longing to the same SExtractor-defined object. In practice 
Gim2D uses thumbnail images, two for each galaxy, ex- 
tracted around the object of interest. The area of these is 
chosen to be 10 times larger than the mask area defined by 
SExtractor. The first thumbnail is cut from the corrected 
frame and is corrected for the local background estimated 
by SExtractor, while the second contains the corresponding 
pixels from the mask image. The fit then produces values 
and uncertainty ranges for each model parameter. Subtrac- 
tion of the best-fit model from the input thumbnail pro- 
duces a residual image that can be used to characterise how 
well the model fits the galaxy. Examples of all these images 



for representative galaxies from each of our morphological 
classes are shown in Figure |5| 

When the fitting algorithm starts to sample the 12- 
dimensional parameter space, it considers not only the pix- 
els assigned to the main galaxy by the mask but all pixels 
flagged as object or background in the SExtractor segmen- 
tation image. Important information about the galaxy could 
be contained in the pixels below the detection threshold. In 
the residual images (see Figure|2J one can see that all pixels 
belonging to other objects in the vicinity of the one under 
consideration are masked out. The final flux is obtained by 
the integration of the best fit model over all pixels, assuring 
that we do not lose the flux in the masked regions. 

The model image of each galaxy is convolved with a 
point spread function (PSF) before comparison with the 
real data. The PSF can be highly variable across a corrected 
frame (StouEhton ct al. 2002) and for this reason it is impor- 
tant to interpolate the PSF parameters measured for indi- 
vidual stars to the position of each galaxy before convolving 
with the model. Our galaxy light model is thus the sum of 
an exponential disk and a Sersic, or de Vaucouleurs, bulge, 
convolved with this "best" PSF. 

Finally a photometric calibration and the redshift are 
required to retrieve physical quantities from the output of 
the B/D decomposition code. 



4 ARE WE CONFIDENT OF OUR 
DECOMPOSITIONS? 

We have carried out the above fitting procedure entirely in- 
dependently for each of our galaxies in the r and the i bands. 
The two segmentation images differ slightly and in addition 
there are colour variations across many of our galaxies as 
a result of variations in the underlying stellar populations 
and in the dust distribution. It is thus reassuring that the 
structural parameters in the two bands are in good agree- 
ment in the great majority of cases (see below). This shows 
that the fitting procedure produces stable results. In addi- 
tion Gim2D produces acceptable converged parameter sets 
for almost all the galaxies in both bands. In Table |21 we 
show the number of galaxies from our spectroscopic sample 
that are successfully modelled by the code in each band in 
the two cases when the Sersic index is set equal to 4 and 
when it is allowed to float. We consider a fit to be success- 
ful when the code is converging. As expected slightly more 
galaxies can be fit when n is kept free, and in this case only 
~ 50 of our 1834 galajcies cannot be fit acceptably in either 
of the two bands. These galaxies are almost all later type 
spirals (Sb-Sc-Sd) which are not modelled simply because 
the centroid position provided by the SDSS database does 
not match the one obtained with SExtractor within the de- 
fined 4 arcsec searching radius. It is interesting to notice that 
the modelling also fails occasionally for early-type objects 
when these are forced to follow a de Vaucolueurs law for the 
central photometric component. We now discuss aspects of 
these fits in more detail. 

4.1 Comparison between fits 

As already noted, we have fitted all our galaxies with a bulge 
model in which the Sersic index n is free and also with a 
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Figure 3. The bulge fraction (upper-left), the disk scalelength 
(lower-left) and the effective radius of the bulge (lower-right) ob- 
tained using different parametric functions to perform the decom- 
position (de Vaucoleurs plus exponential on the x-axis and Sersic 
plus expenential on the y-axis) are plotted against each other for 
1702 galaxies in the i band. 

model in which it is fixed to the de Vaucouleurs value n = 
4. In this subsection we show that the parameters of the 
decomposition of most interest to us are only weakly affected 
by this choice for most galaxies. In Figure|3|we compare the 
values of the bulge fraction B/T, of the disk scale length h 
and of the bulge effective radius re obtained for each galaxy 
in the two cases. We show results for the r band only (results 
for i are similar). 

4.2 Comparison between different bands 

In this section we compare the parameters estimated for each 
galaxy when the same model is fit independently to images 
in each of the two SDSS bands analysed here. In Figures 2] 
and|S]we plot the retrieved parameters in the Sloan i and r 
bands against each other for the 1636 galaxies modelled suc- 
cessfully in both bands assuming a de Vaucouleurs profile for 
the bulge and an exponential for the disk. The correspond- 
ing results are shown in Figure |S| and |7| for the 1766 galaxies 
modelled successfully in both bands assuming a Sersic pro- 
file for the bulge and an exponential for the disk; in this 
case, however, we also compare the values of the index n 
found in the two cases. 

In the apparent magnitude plot there are only few iso- 
lated points while the rest correlate very well, the scatter 
being consistent with the expected variation in mean colour. 
The scatter in bulge-to-total light ratio is also gratifyingly 
small for the vast majority of the galaxies, particularly when 
n is fixed to the de Vaucouleurs value. The scatter is slightly 




Figure 4. The total flux (upper— left), the bulge fraction (upper- 
right), the disk scalelength (lower— left) and the effective radius of 
the bulge (lower-right) in the i and r bands are plotted against 
each other for the 1638 galaxies modelled in both bands by fitting 
a de Vaucouleurs profile to the bulge and an exponential to the 
disk. 

larger when n fioats because difi'erent best fit values of the 
index in the two bands (see Figure |7J lead to different splits 
of the luminosity between bulge and disk. In most cases, 
however, similar n values are found in the two bands. The 
"arrow" shape of the B/T plots reflects the fact that the 
data in one band occasionally prefer a weak disk while no 
disk is present in the best fit in the other band. 

There is also quite good agreement between the values 
of the scalelengths for the bulge and disk components mea- 
sured in the two bands. The agreement is worse for bulges 
than for disks and for components of small angular size com- 
pared to larger ones. This is presumably a refiection of reso- 
lution problems due to the finite pixel size and to difficulties 
with the PSF deconvolution. Nevertheless the apparent axial 
ratios of both bulges and disks agree well in the two bands 
with the scatter increasing for rounder systems. The agree- 
ment between the position angle of the disk in the considered 
photometric bands is pretty good even if a substantial scat- 
ter is present. The points in the two extreme corners could 
be moved by a simple rotation of 180° . In addition we check 
that the distribution of the position angle is consistent with 
a random orientation on the sky. 

4.3 Error estimates and goodness of fit 

When carrying out its fitting Gim2D constructs a value 
for each PSF-convolved model by summing over all pixels 
within the mask the square of the difference between model 
and data divided by the variance of the pixel noise, assumed 
to be due entirely to photon statistics. This measure of good- 
ness of fit is then minimised over all parameters (each re- 
quired to lie within a prespecified "allowed" range) to lo- 
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Figure 5. The position angle of the disk (upper— left), the ellip- 
ticity of the bulge (lower— left) and the inclination angle of the 
disk (lower— right) in the i and r bands are plotted against each 
other for the 1638 galaxies modelled in both bands by fitting a de 
Vaucouleurs profile to the bulge and an exponential to the disk. 




Figure 6. The same as Figure 111 but for the 1766 galaxies mod- 
elled in both bands by fitting a Sersic profile to the bulge and an 
exponential to the disk. 

cate the maximum likelihood model. Once the algorithm 
has converged, the region of parameter space surrounding 
the likelihood maximum is sampled in order to compute 
marginalised a posteriori one- dimensional probability dis- 
tributions for each model parameter. These are then used 
to define best parameter estimates, taken to be the medians 



Figure 7. Similar to Figure|21but for the 1766 galaxies modelled 
in both bands by fitting a Sersic profile to the bulge and an ex- 
ponential to the disk. In this figure we also compare the values of 
the Sersic index n found in the two bands (upper right). 



of these distributions, and 99% confidence ranges defined 
by their upper and lower 0.5% points. Nevertheless he com- 
puted turns outs to be not too sensitive to whether prob- 
lems occur in the decomposition (i.e. a wrong point spread 
function) or the decomposition reliably describes the hght 
distribution in the galaxy. 

With the aim of better understanding the goodness of 
the fit of our models we introduce two additional measures, 
Gi and G2, which characterise the size of the residuals with- 
out reference either to the overall luminosity and size scales 
of a galaxy or to the expected counting noise. They are de- 
fined using the region flagged as belonging to the galaxy 
in the segmentation image generated by SExtractor. Gi is 
the difference between the model and observed luminosities 
in this region as a fraction of the model luminosity, while 
G2 is the ratio of the sum of the absolute values of the 
residuals in all the pixels to the model luminosity. Work- 
ing from the individual pixels (ij), by definition the total 
counts in the science image Oij) and in the model image 
Mij) are due to the hght of the galaxy plus a uniform 
sky (Y^ Sij), which is the same in the two images. The total 
counts in the residual image simply reflect the difference be- 
tween the luminosity of the observed and modelled galaxy 
(Y^ Dij = Oij—Y^ Mij)- Our deflnitions can consequently 
be formulated as: 

Gi = ^^ and G2 = ^^ (3) 

Figure |H| shows the distribution of galaxies with respect to 
these quantities as a function of galaxy absolute magnitude. 
Results are presented for flts to the i band data both for 
floating n and for n fixed to the de Vaucouleurs value. Re- 
sults for the r band are very similar. The Gi parameter is 
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Figure 8. Goodness of fit measures for our two-component mod- 
eUing of the j-band images. The two left-hand panels give results 
for the 1469 galaxies modelled successfully using a de Vaucouleurs 
profile for the bulge and an exponential profile for the disk. The 
two right-hand panels give results for the 1532 galaxies modelled 
successfully when we use the more general Sersic profile for the 
bulge. The two upper panels show results for Gi , the total light in 
the residual image in units of the total light of the model. The two 
lower panels give results for G2 , the sum of the absolute values of 
the individual pixel deviations from the model again in units of 
the total model luminosity. Gi and G2 are plotted against abso- 
lute magnitude in the main part of each panel while a histogram 
of their marginal distribution is given in the inset. 



narrowly distributed around zero, with a slight bias towards 
positive values. Thus Gim2D underestimates slightly the 
luminosities of these large galaxies, but typically by only a 
couple of percent. Misestimates by more than 5% are very 
rare. The distribution of G2 peaks at 0.1 and is skew with a 
longer tail towards higher values. There is a tendency for de- 
viations from the models to be larger for intrinsically fainter 
galaxies, particularly below about Ali = —20. Residuals are 
only slightly reduced by the extra freedom involved in al- 
lowing n to vary because, as can be seen from the examples 
in Figure ^ the dominant residuals are often due to non- 
symmetric features such as spiral arms or dust lanes. We 
are encouraged that G2 is less than 15% for about three 
quarters of our galaxies and almost never rises as high as 
25%. We therefore believe that our model represents the im- 
ages of the majority of galaxies adequately for our purpose, 
and that derived parameters can be used meaningfully to 
characterise physical properties of the galaxies themselves. 
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Figure 9. Distribution of the Scrsic index n. The upper-left panel 
shows the distribution for the total sample of 1790 and 1782 galax- 
ies successfully modelled in the i and r band respectively. In the 
other panels the Sersic index distribution is separated according 
to morphological type. 



well. While massive systems usually require large values of 
n, similar to the de Vaucouleurs value, less massive ellipti- 
cals and bulges, particularly the bulges of late-type spirals 
usually demand smaller n values and indeed can often be 
well fit by an exponential law with n — 1. For our JPG 
sample when we allow the Sersic index to assume any value 
between 0.2 and 4 we obtain the distributions shown in Fig- 
ure |5| The upper-left histograms give results for the 1790 
galaxies in the i band and the 1782 galaxies in the r band 
(out of the the total sample of 1862) which are successfully 
fit by Gim2D. These distributions peak at n = 4, confirm- 
ing that the r^^'^ law provides an acceptable fit to the bulge 
component of a large fraction of the galaxies in a magnitude- 
limited sample. In the other panels of Figure |5| we split the 
sample by morphological type and it becomes evident that 
most E/SO galaxies are well described by values of n close to 
than 4. The same is true for the bulges of most early-type 
spirals. For later-type spirals there is a clear shift to lower 
values of n. for many of these objects n 1 is preferred, 
confirming the earlier studies referred to above. Note that 
the distributions in Figure |5] are almost independent of the 
band in which the decomposition is carried out, confirming 
the robustness of the results. 



5 RESULTS 

5.1 Comparison of Sersic index distributions 

Previous studies jAndredakij Il998l : ICourteau et"ai] Il996l: 
Ide JonJl99^) have shown that the Sersic formula can fit the 
light profiles of many nearby ellipticals and bulges extremely 



5.2 Disk and Bulge luminosity 

In this section we derive the fraction of the luminosity den- 
sity in the local universe in bulges and in disks for the Sloan i 
and r bands, starting from our complete sample of 1 — selected 
galaxies with r < 15.9 (after correction for Galactic extinc- 
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Galaxies modelled 


Type of fit 


r band 


i band 


de Vaucoulcurs + exponential 
Sersic + exponential 


1450 
1528 


1469 
1532 



Table 2. Galaxies in the spectroscopic sample that we use to calculate the fraction of light in the local universe in disks and bulges. Of 
the 1588 galaxies in our sample for which the redshift is known, 1517 are successfully modelled by the code in both bands when using 
the Sersic profile and 1409 if we adopt the de Vaucouleurs law. In the table we show the number of galaxies in the i and r bands for 
which we have the spectroscopy, the decomposition parameters and good photometry. 



Typo of fit 


r band 


i band 


de Vaucouleurs + exponential 
Sersic + exponential 


(58.92 ±2.40)% 
(54.82 ± 1.95)% 


(54.92 ± 2.02)% 
(55.41 ± 1.98)% 



Table 3. Total fraction of the light in disks in the local universe in the r and i bands and for fits requiring n = 4 and allowing n to vary 
over 0.2 < n < 4. 
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Figure 10. Distribution of apparent axial ratio for the disk com- 
ponents of the 1469 galaxies with a redshift which were modelled 
successfully in the i band by a de Vaucouleurs plus an exponential. 
Galaxies are split by absolute magnitude into 10 bins containing 
approximately equal numbers of objects. 



tion). From our total sample of 1862 galaxies we here con- 
sider only the 1588 objects for which spectroscopic data are 
available and it is therefore possible to measure the abso- 
lute magnitudes needed by our estimation procedure. The 
absolute magnit udes used in this pap er are k-corrected us- 
ing the code of lBlanton et alJ l|20ol), i;2_16. In the SDSS 
main galaxy sample as a whole the median redshift is near 
z = 0.1, so Blanton chose to express results in the SDSS 
filter system shifted by 0.1. For our sample the median red- 
shift is about 0.05. Nevertheless for consistency with other 
SDSS work (in particular, with the luminosity functions we 
use below) we follow Blanton's convention and k-correct to 
z — 0.1. We denote absolute magnitudes in this system as 



Figure 11. As Figure ITol but showing the distribution over mor- 
phological type. 

"'^Mr and "'^M; to distinguish them those in the unshifted 
system. Since our primary results concern the ratios of lu- 
minosities in different components, this choice has no effect 
on our analysis 

There are 23 objects in this set for which no morpho- 
logical type was assigned by the JPG astronomers. We ex- 
clude these from further consideration here, leaving 1565 
objects with redshift and a well defined "by eye" morpho- 
logical type. From this sample we also excluded all galaxies 
for which Gim2D failed in the modelling, and in addition we 
removed three objects for which the k-correction is not reli- 
able due to bad photometric data in the bluest and reddest 
bands. This reduced the number of objects used to estimate 
the luminosity densities in bulges and disks to the numbers 
in Table H 

Our strategy for computing the fraction of the local lu- 
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Type of fit 


r band 


i band 


de Vaucouleurs + exponential 
Sersic + exponential 


(33.40 ±0.75)% 
(34.39 ±0.72)% 


(31.15 ±0.73)% 
(35.82 ±0.74)% 



Table 4. Total fraction of the light in pure bulge systems in the local universe in the r and i bands and for fits requiring n = 4 and 
allowing n to vary over 0.2 < n < 4. 
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Figure 12. As Figure lTni but showing the distribution over bulge- 
to-disk ratio. 



minosity density which is in bulges and disks is as follows. 
We separate the galaxies in our sample into 11 bins accord- 
ing to their absolute luminosities in the r or i bands. The 
nine brightest bins each contain about 10% of the sample 
while the two faintest bins contain about 5% (we made this 
choice in order to get better luminosity coverage for faint 
galaxies.) For each bin we then use our decompositions to 
estimate the fraction of the total light coming from disks for 
galaxies at that absolute magnitude. Assuming this fraction 
to be appropriate for all galaxies of similar intrinsic bright- 
ness, we combine it with luminosity functions determined 
from much larger SDSS samples by Blanton ct al. (2003) to 
obtain the fraction of the local luminosity density which is 
in disks. The complementary fraction is then the amount in 
bulges. 

While this appears straightforward, a serious complica- 
tion arises from the fact that 2-D fitting codes like Gim2D 
tend to fit radial variations in axial ratio or position angle 
in ellipsoidal galaxies by assigning a fraction of their light 
to a disk, when in fact none is prese nt. This systematic i s 
well known and is commented on in ISimard et all J200d) . 
We can demonstrate its presence in our purely luminosity- 
selected data by examining the distribution of disk apparent 
axial ratio returned by Gim2D. The distribution of hja is 
expected to be uniform on [0, 1] for randomly oriented thin 
disks. Figure ITUl shows the distributions we actually obtain 



for the "disks" in our sample, split into 10 equal bins by 
absolute total i-magnitude. While for faint galaxies these 
distributions are indeed consistent with being flat, in the 
brighter bins there is clearly a strong bias towards high b/a. 
Among the brightest galaxies almost no "disks" are found 
with hja < 0.5. Figure fTTI shows that these bright bins are 
dominated by early-type galaxies according to the visual 
classifications of the JPG. The absence of small b/a values 
demonstrates that few of these systems actually have signif- 
icant thin disks, despite the fact that Gim2D assigns most 
of them B/T ratios substantially smaller than unity (see 
Figure IT^. 

In order to circumvent this problem in the following 
analysis we use only the galaxies in each absolute magnitude 
bin which have b/a < 0.5. We believe the great majority of 
these must be true disks since ellipticals with apparent axial 
ratios smaller than 0.5 are very rare. For random orientation 
the total number of true disks expected in the bin is just 
twice the number with b/a < 0.5. The total light in disks in 
the bin is, however, more than twice the light in the disks 
with b/a < 0.5, since dust extinction is significantly stronger 
in edge-on than in face-on disks. This must be corrected if 
we wish to obtain an unbiased estimate of the luminosity 
density in disks. 

In practice, our procedure works as follows. For each 
absolute magnitude bin k we estimate the fraction of the 
light in the disk component as 



fdisk,k 



Ldis 



Ltot.k 



(4) 



where Ldisk,k is the total luminosity of the disks of the galax- 
ies in the bin and Ltot,k is the total luminosity from all com- 
ponents of these same galaxies. Since we assume that we can 
rely on our decomposition only for edge-on systems, we split 
the numerator into two parts 



.fdu 



L. 



b/ a<^0.5,disk,k 



+ L, 



b/ a>Q.5,disk,k 



Ltot.k 



(5) 



where Li,/a<o.5,diBk,k is the luminosity due to "edge-on" 
disks with b/a < 0.5 and is obtained directly from our de- 
compositions. We estimate the luminosity Lf,/a>o.5,disk,k in 
"face-on" disks by assuming that true disks are randomly 
oriented and that their internal extinction A\ depends on 
inclination according to the standard prescription 



Ax = 7Alog(a/6). 



(6) 



Here Ax is the correction to exactly face-on orientation. 
Assuming this formula, it is a simple matter to relate the 
total luminosity density in disks to that in disks with b/a < 
0.5. We find 



Ldisk.k — 2 



1 + 0.47;^ 



hence 



J disk, k — ^ 



1 + 0.47;^ -Lb/a<0.5,disk,k 



b/a<Q.5,disk,k^ (7) 

(8) 



Ltot,k 
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The bulge to total ratio B/T is one of the structural param- 
eters returned by Gim2D for each galaxy. In combination 
with the total absolute Petrosian magnitude (taken directly 
from the SDSS database) it allows us to estimate the disk 
luminosity of each object contributing to ii,/a<o.5,iiisfc,fc- We 
obtain Ltot,k simply by summing the Petrosian luminosities 
of all galaxies in the bin regardless of their b/a. We take 
values for 7a from the work of TuUy et al. (1998). Since the 
numerical coefRcient on the right hand side of equation |H] 
is only slightly smaller at i than at r band, we assume the 
value 2.56 for both photometric bands. 

The results of these calculations are shown in Figure ITHl 
for the galaxies modelled successfully in the i and r bands. 
In each figure we show results separately for decompositions 
which force a de Vaucouleurs bulge and for decompositions 
in which we allow 0.2 < n < 4. The pattern is very similar in 
all cases. The average light fraction in disks varies smoothly 
from about 10% in the brightest galaxies to almost 100% 
in faint galaxies. For the brighter bins these fractions are 
much smaller than the values obtained from a direct naive 
average of the B/T histograms of Figure m\ because of the 
systematic effect we have just been discussing. 

The error bars on the points in Figure [T^ are important 
because they determine the precision of our final results. 
We neglect the formal errors on B/T returned by Gim2D 
because these are quite small, typically ±7%, and are well 
below the systematic error discussed above due to the as- 
signment of isophote twists or axial ratio changes to spu- 
rious thin disks. We assume, however, that this systematic 
can be neglected for systems with b/a < 0.5. The uncer- 
tainty in our estimate of the disk light fraction in each bin 
is then dominated by sampling. As a statistical model for the 
population of a particular bin we assume that a randomly 
chosen galaxy has a detectable edge-on disk {b/a < 0.5) 
with probability p^o where a priori we have < Peo < 0.5. 
We also assume that the B/T values of these edge-on sys- 
tems are drawn at random from some unknown distribution 
with population mean and variance which we estimate us- 
ing the sample mean and variance, {B/T) and Var(_B/r) re- 
spectively. If the bin contains Nt galaxies of which N^o have 
disks with b/a < 0.5, then the maximum likelihood estimate 
of Peo is Peo = Neo/Nt provided Neo/Nt < 0.5 (it is equal to 
0.5 otherwise). To approximate the variance of peo we use 
the standard binomial formula Var(peo) — Peo(l — Peo) / Nt, 
even though this is formally incorrect for peo ~ 0.5. Our 
estimate of the mean light fraction in edge-on disks is then 
(1 — {B/T))peo and we calculate the variance in this esti- 
mate as {peofVw{B/T) + {B/TfYar{peo)- These give the 
final results plotted when multiplied by the correction factor 
of equation |7| which accounts for the light in face-on disks. 
There is undoubtedly a systematic uncertainty associated 
with this last step, but we ignore it here. 

We can now average the disk light fractions of Figure 
1131 over the galaxy population as whole in order to obtain 
the fractions of the total luminosity density in the local uni- 
verse coming from disks and from bulges. The contribution 
of each of our absolute magnitude bins to the total luminos- 
ity de nsity ^tot.k c an be obtained by integrating the appro- 
priate |Blanto^.etjil. (2003) luminosity function across the 
bin. The final result for the fraction of the local luminosity 



density in disks is then, 

J, y ] ^ ^tot,fc ' fdisk,k 

Jdisk — ^ • (,9j 

2^4, ^tot,k 

Note that the luminosity fun ctions for our sample are quite 
consistent with those given bv lBlanton et all i2003fl . We pre- 
fer to use the latter here because of their much better sta- 
tistical precision. The final result we obtain for the total 
fraction of the light coming from disks in the local universe 
is (54±2)%, with no detected dependence on observing band 
or decomposition parametrisation. The details are in Table 
13 The error bars in this table are calculated directly from 
those in Figure lT^ assuming the uncertainties in the different 
absolute magnitude bins to be independent. Uncertainties 
coming from the luminosity function itself are negligible in 
comparison. 

A slight variation of this analysis allows us to calculate 
a second interesting quantity: the fraction of galaxies in each 
of our absolute magnitude bins which contain no detectable 
thin disk and so may be considered "pure" bulge systems. 
Our hypothesis here is that Gim2D will detect any signifi- 
cant disk if it is sufficiently inclined to the line-of-sight that 
b/a < 0.5. Exactly one half of all disk galaxies should be at 
least this inclined. We can thus estimate the number of effec- 
tively diskless galaxies in each bin by subtracting twice the 
number of objects with b/a estimates below 0.5 from the to- 
tal number of galaxies in the bin. The result of this exercise 
is shown in Figure ITU in similar format to Figure IT^ Again 
the results are very similar in the two pass-bands and for our 
two assumptions about bulge profiles. For the brightest bins 
we find that the great majority of galaxies are effectively 
diskless, while in the faintest bins our statistics are consis- 
tent with at most a small fraction of "pure bulge" systems. 
The fraction of "pure bulges" varies smoothly with absolute 
magnitude between these two extremes. The error bars re- 
flect the binomial uncertainty in our estimate Peo ~ Neo/Nt 
of the fraction of the bin population with detectable edge-on 
disks. 

As before we can combine the results of Figure IT^ with 
the luminosity functions of Blanton et al. (2003) to estimate 
the fraction of the luminosity density of the local Universe 
which is contributed by effectively diskless systems. The re- 
sult, given in detail in Table |1| is about (32 ± 1)% with no 
significant dependence on pass-band or bulge fitting func- 
tion. Thus the breakdown of stellar luminosity in the local 
Universe is apparently 54% in disks, 32% in "pure bulge" 
systems with no photometrically detectable (by Gim2d) 
disk and only 14% in the bulges of galaxies with detectable 
disks. 



6 CONCLUSIONS 

We have used the tw o-dimensional p hotom etric fitting pro- 
gramme Gim2D of ISimard et alJ i2002f) on the images 
of an apparent magnitude limited sample of 1862 large 
galaxies selected from the SDSS and visually classified by 
Jj^akamura et al. (2003^. In almost all cases the code returns 
a well-defined decomposition of the galaxy into disk and 
bulge components with parameters which have small formal 
error bars and vary little either between the r and i band 
images analysed here or between decompositions in which 
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Figure 13. The two left-hand panels show the luminosity frac- 
tion in disks for galaxies in the Sloan i band. The galaxies are 
modelled with two different sets of parametric functions: a de 
Vaucoulours profile for the bulge and an exponential for the disk 
(left panel) or a Sersic profile for the bulge plus an exponential 
for the disk (right panel). The two right-hand panels show the 
same results but in the r band. 



Figure 14. The two left-hand panels show the fraction of pure 
bulge galaxies as a function of the i band absolute magnitude. 
Different parametric functions are considered: de Vaucouleurs 
plus exponential (left— panel) and Sersic plus exponential (right- 
panel). The two right-hand panels show the same results but in 
the r band. 



the Sersic index of the bulge component varies or is fixed 
to the de Vaucouleurs value n — 4. The total amount of 
light in differences between the observed image and the best 
fit model is typically only about 10% of the galaxy lumi- 
nosity. Despite this apparent success, we show that for in- 
trinsically bright galaxies the "disk" component does not in 
most cases represent a true thin axisymmetric disk, since the 
sample distribution of axial ratios b/a deviates strongly from 
the uniform distribution expected for a population of such 
disks - near edge-on disks are g rossly under rcprcscnted. This 
problem was already noted bv lSimard et al. (2002,) . Appar- 
ently the code is using the degrees of freedom provided by 
the assumed disk component to fit radial changes in isophote 
shape or position angle within galaxies which have no real 
thin disk. 

We attempt to correct for this systematic problem by 
concentrating on galaxies for which Gim2D finds a disk 
axial ratio b/a < 0.5. We argue that ellipsoidal stellar sys- 
tems with such extreme apparent axial ratios are quite rare 
so that the component isolated by the the decomposition 
programme is likely to be a highly inclined disk. Since the 
selection criteria for our sample depend at most weakly on 
inclination, we can use the assumption that our galaxy sam- 
ple is randomly oriented to correct from the sample with 
highly inclined disks to the sample as a whole. Such correc- 
tions are at best approximate, since there are undoubtedly 
cases where the decomposition mixes hght from the true 
bulge and disk components in such a way that photomet- 



rically fit "disk" has b/a > 0.5 even though the galaxy is 
sufficiently inclined for the projected image of its true disk 
to have b/a < 0.5. Without further information (for exam- 
ple, from kinematics) it is very difficult to assess the size of 
the biases this may introduce. 

As a first application we estimate average disk light frac- 
tions for galaxies as a function of absolute magnitude. These 
range from about 10% for the brightest galaxies to almost 
100% for the faintest ones. At each absolute magnitude we 
also estimate the fraction of "pure bulge" galaxies, defined 
as galaxies for which Gim2D would detect no disk with 
b/a < 0.5 even if the orientation were such that b/a < 0.5 
would be expected for any true axisymmetric thin disk. We 
find that most of the galaxies in the brighter absolute mag- 
nitude bins are "pure bulge" by this definition, but that 
this fraction decreases steadily for fainter systems. We do 
not detect any population of "pure bulges" in our faintest 
bin. These numbers differ substantially from those inferred 
naively from the B /T distributions measured directly by 
Gim2D: for example, in the brightest absolute magnitude 
bin about 85% of galaxies are "pure bulge" and disks only 
contribute ~ 10% of the total Ught, yet the median value of 
B/T returned by the code is 0.55 and very few galaxies are 
fitted with B/T > 0.8 

By combining these results for the absolute magni- 
tude dependence of the mean light fractions in disks or in 
"pure bulges" with the luminosity functions of lBlanton et alJ 
i2003r . we have been able to estimate the fractions of all 
galaxy light in the local Universe coming from disks and 
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from "pure bulges". The results depend little on the wave- 
band used or on the bulge luminosity profile we adopt. We 
find that 54 ± 2% of the local luminosity density is con- 
tributed by stars in disks and 32 ± 2% by stars in "pure 
bulges" . The remaining f 4 ± 3% comes from bulges in sys- 
tems with detectable disks. The mean bulge-to-total ratio of 
the latter systems is thus 14/68 ~ 0.2, substantially smaller 
than typical B /T values in the histograms of FigurelH Note 
that only sampling uncertainties are taken into account in 
the errors quoted here. Residual systematics may remain, 
reflecting the fact that real galaxies have a more complex 
structure than the models we use to describe them here. A 
fully convincing decomposition of galaxies into bulge and 
disk components clearly cannot be performed using imaging 
data alone. We believe, however, that the analysis of this 
paper provides the best quantitative indication so far of the 
overall distribution of stars between the two basic structural 
forms which make up galaxies. 
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